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Q> ■ ABSTRACT 

We construct spherical, hydrostatic models of dense molecular cores and Bok globules consisting of two dis- 
tinct, spatially separate gas components: a central, isothermal region surrounded by a negative-index, polytropic 
bJ)| envelope. The clouds are supported against their own self-gravity by a combination of thermal, mean magnetic, 

^ i and turbulent wave pressure. The latter two are included by allowing for locally adiabatic, non-isentropic pressure 

' components. Such models are meant to represent, in a schematic manner, the velocity and density structure of 

cores and globules, as inferred from molecular line and dust continuum observations. In addition, our picture 
reflects the theoretical expectation that MHD wave motions, which are important at scales ^0.1 pc in typical 
low-mass star-forming regions, are damped at smaller scales, giving rise to a finite-sized, thermally-dominated 
core region. We show that if the pressure components are isentropic, then the pressure drop from the center to 
the edge of the composite polytropes we consider is limited to 197, the square of the value for the Bonnor-Ebert 
I sphere. If the pressure components are non-isentropic, it is possible to have arbitrarily large pressure drops, in 

. agreement with the results of McKee & Holliman (1999). However, we find that even for non-isentropic pressure 

CO ' components, the ratio of the mean to surface pressure in the composite polytropes we consider is less than 4. We 

, show by explicit construction that it is possible to have dense cores comparable to the Jeans mass embedded in 

stable clouds of much larger mass. In a subsequent paper, we show that composite polytropes on the verge of 
gravitational instability can reproduce the observed velocity and density structure of cores and globules under a 
variety of physical conditions. 

Subject headings: ISM: clouds — ISM: globules — ISM: structure — Stars: formation 
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|H , L INTRODUCTION thermal motions are responsible for at least half of the observed 

c/2 ■ „, .. f iJ J 11 1 J ■ • . line broadening. This is in contrast to the weU-known increase 

rri . Observations of cold, dense molecular clouds in the inter- ^ , ,°. , . , 

. „ . . , „ » »u • » T 1 » J of nonthermal linewidth with size scale in the regions imme- 

• • ' stellar medium present a unique challenge to theorists. Isolated ... ,. , . n 

D,,,, 1, . iiij diately surrounding the ammonia cores (Fuller & Myers 1992; 

Bok globules and dense molecular cores in molecular cloud _ / , ^^Jl, _, . ^ . ^ 

1 . ,.1 J £ . J f • f Goodman et al 1998). The properties of cores in more crowded 

complexes, in particular, are definitive testbeds for theories of , , ^ ,. , . . ■ r„ 

111,,, . 1 , f . regions, such as Ophiuchus and Orion, differ from those in Tau- 

« molecular cloud structure and low-mass star formation. As , , 

Z3 I ,,,, . 1 , , , • J • . . 11 J- " rus mainly insofar as they tend to be more massive and compact, 

the simplest unit structures in the dense interstellar medium , , . , -' , ,. . , , ,^ „• n . , 

■ ■ .j^- 1 o r^i inoT\ J .J • f and have higher temperatures and linewidths (Casein & Myers 

(Dickman & Clemens 1983), an adequate description of these . , , n ^, ■ i^^^n 

\ . , . . u . ■ .1. X . A- f 1995; Motte, Andre, & Neri 1998). 

objects should contribute significantly to our understanding of ^ , ^ / 



■ 
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^,11 , 1 u 1 A catalog of 248 optically selected Bok globules was com- 

the molecular component of the Galaxy as a whole. ., ,, _ J; ,,„„„, . ° 

piled by Clemens & Barvaims (1988). A comparison of these 

objects with the Benson & Myers (1989) sample of dense cores 

1.1. Observed Properties of Dense Cores and Bok Globules molecular cloud complexes reveals that the two populations 

Dense cores are compact regions within molecular clouds are quite similar, with a trend toward higher masses and den- 

whose number density exceeds about 10"* cm'l The connec- sities (by about a factor of two) in the dense cores (Bourke et 

tion of these regions to star formation in the complexes has al. 1995). Associated IRAS point sources have been detected 

been explored in some detail (Loren 1989; Onishi et al 1998); i" about one-quarter of the Bok globules, down to a detection 

in many cases there is close spatial agreement between young limit of 0.7 Mq (Yun & Clemens 1990). The degree of sim- 

stars and cores (Beichman et al 1986). In fact, IRAS sources ilai'ity is perhaps surprising, given the different environments 

have been detected in approximately half of the known (~ 100) in which the two types of object reside. For example, the two 

nearby cores (Benson & Myers 1989). The relatively isolated populations should differ with respect to their exposure to the 

cores in Taurus, which have been extensively surveyed in NH3 interstellar radiation field. Whereas cores in complexes are of- 

(Myers & Benson 1983; Benson & Myers 1989), have sizes ten shielded from such radiation by their parent clouds (Ben- 

on the order of 0.1 pc, masses - 1 - 10 Mq, kinetic tempera- son & Myers 1989; Onishi et al 1996), photoelectric heating 

hires ~ 10 K, and NH, linewidths ~ 0.3 km s"'. The narrow ^om the intersteUar radiation field is expected to produce an 

Unewidflis suggest that these objects are relatively quiescent; ouWaid temperature gradient in Bok globules lacking embed- 

'Present address: Department of Physics and Astronomy, University of Western Ontario, London, ON N6A 3K7 
^Email; curry@astro.uwo.ca 
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ded young stars (Boland & de Jong 1984). This effect has been 
observed in dark clouds and cores in close proximity to star- 
forming regions (Young et al 1982; Bachiller, Guilloteau, & 
Kahane 1987), but direct evidence for isolated Bok globules is 
more scarce. According to Clemens, Yun, & Heyer (1991), the 
spectral energy distributions of globules in the Clemens & Bar- 
vainis sample preclude single-temperature models for all ob- 
jects. On the other hand, molecular line observations of objects 
in the same sample display remarkably uniform temperatures 
over their projected areas (Lemme et al 1996). The latter might 
be explained, however, by the fact that the dense molecular trac- 
ers used in such studies do not sample the outer, more exposed 
regions of the globules. 

1.2. Virial Parameters: Implications for Stability 

A key question relating to cores and Bok globules is: how 
near are they to gravitational collapse? A stable interstel- 
lar cloud need not be highly centrally condensed and self- 
gravitating; it can also be nearly uniform and confined by an 
external pressure. In an analysis of clump data ^ from several 
active star-forming regions, Bertoldi & McKee (1992) noted 
a correlation between the "observed virial parameter" for a 
clump, a = 5R{a^) /GM, and the clump mass, M. Here, (ct^) 
is the mean square velocity along the line of sight and R an av- 
erage cloud radius, both determined from intensity contours in 
the 7 = 1 - Hne of ^^CO."^ Stable clumps with a > 1, they ar- 
gued, are likely to be confined by the pressure of the external 
medium, and not very centrally concentrated, while those with 
a ^ \ would be expected to be strongly condensed and self- 
gravitating (no clumps with a<^\ were observed; such clumps 
are unstable and would collapse unless supported by a strong 
magnetic field). The authors' analysis revealed that only the 
most massive clumps (;> 100 Mq) in the complexes had a ^ 1, 
and thus that the majority of the clumps seen were likely to 
be pressure-confined. For example, representative values for a 
^ IOMq clump in Ophiuchus from Loren (1989) are ~ 0.46 
km s"' and R ~ 0.43 pc, implying a ~ 10.5. 

It is important to emphasize that estimates of a from molecu- 
lar line studies are highly dependent on the observed molecule. 
Dense molecular tracers, such as NH3 and CS, are generally 
only excited in small regions where the mean density of hydro- 
gen exceeds ~ 10"* cm"^. The mean density of the Ophiuchus 
'"^CO clumps, on the other hand, is only ~ 3 x 10^ cm"^. It is 
therefore desirable to recalculate a for these denser molecular 
tracers, in order to ensure that one is truly probing the densest 
gas. For this purpose, we turn to the NH3 study of dense cores 
by Myers & Benson (1983). Mean values from their sample 
are: a ~ 0.23 km s"', ^ ~ 0.06 pc, and M ~ 4.1 Mq. These 
yield a ~ 0.9, indicating that the cores in their sample, which 
are of very low mass compared to those in the '^CO studies, 
are nevertheless nearer to a state of marginal stability. Further 
support is provided by dust continuum studies in Ophiuchus, 
which also give a ^ 1 for cores of this size and smaller (Motte 
et al 1998). 

In this paper, we adopt the view that both low- and high- 
mass cores that have been detected in ammonia and other 



dense molecular tracers (e.g. C'**0 and CS) are near a state of 
marginal stability, or, as we shall refer to it, near the critical 
point. The latter is more precisely defined as the unique state 
such that, upon suffering a small perturbation at its surface, the 
cloud must collapse. These cores may themselves be embedded 
in molecular gas that is gravitationally bound to the core; how- 
ever, we assume that this external gas acts primarily to set the 
pressure at the surface of the core, not to determine the cloud 
stability. 

1.3. The Initial Conditions for Star Formation 
1.3.1. Pre-Protostellar Cores 

The high percentage of dense cores and Bok globules with 
embedded IRAS point sources can, in principle, put restric- 
tions on theories of star formation. Furthermore, the roughly 
equal samples of "starless" cores and Bok globules should place 
constraints on the initial conditions for protostellar collapse, 
assuming that these clouds still have a reasonable chance of 
forming stars (§1.2). Hence the recent usage of the term "pre- 
protostellar core" to denote these objects (Ward-Thompson 
et al 1994; Andre, Ward-Thompson, & Motte 1996). We 
shall henceforth refer to both starless dense cores found within 
GMCs and starless Bok globules by the latter term, for brevity 
abbreviated to PPC. In the quasi-static picture of low-mass 
star formation, PPCs evolve from low density and possibly 
low mass along a near-equilibrium sequence of denser, pos- 
sibly more massive states, until a state of marginal stability 
is reached. This evolution may be achieved in a number of 
ways: ambipolar diffusion is likely to play a key role, and pos- 
sibly ongoing accretion from the surrounding medium (Mestel 
& Spitzer 1956; Boland & deJong 1984; Bonnell, Bate, & Price 
1996). As defined, PPCs generally only form one or at most a 
small group of stars. 

Historically, theories of cloud equilibria have been intimately 
connected with calculations of gravitational collapse, inasmuch 
as the former fill the important role of initial conditions for the 
latter. Models with singular density profiles are often adopted 
for computational simplicity. In particular, the singular isother- 
mal sphere (hereafter SIS), which has a static density structure 
given by p{r) cx r~^, has been extensively studied. This model 
is unstable to collapse, and has a corresponding time-dependent 
solution first studied by Shu (1977). 

As a model for PPCs, there are three principal caveats to the 
SIS, two of which stem directly from observation. The first 
is the existence of a nonthermal component of the molecular 
linewidth in PPCs. The nonthermal line width increases with 
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size scale as r in low-mass cores and as r in hig h-mass 
cores (Caselli & Myers 1995). 
thermal motions prevail within 
imum 



As mentioned in §1.1 



0.1 pc of the intensity max- 
but even along the line of sight toward the intensity 
maximum, a finite nonthermal component is detected (Good- 
man et al 1998). Second, most objects that have been mapped 
with sufficient resolution exhibit non-singular — indeed, nearly 
flat — column density profiles as one approaches the intensity 
maximum. This is true of both stable Bok globules (Dickman 
& Clemens 1983; Turner, Xu, & Rickard 1992; Lehtinen et 



'We shall follow the terminology of Williams, Blitz, & McKee (1999) in describing the structure of molecular clouds. "Clumps" are coherent regions in l-b — v 
space, generally identified from spectral line maps of molecular emission. Star-forming clumps are the massive clumps out of which stellar clusters form; cores 
(sometimes refen'ed to as "dense cores") are the regions out of which individual stars or stellar systems form. The clumps studied by Bertoldi & McKee were iden- 
tified based on '^CO line maps, whereas dense cores in low-mass star-forming regions are generally within star-forming clumps and are identified from their NH3 
emission. 

''Assigning a "radius" to a non-spherical clump or core is common practice. The authors showed that, if a is interpreted as as average over all orientations of an 
ellipsoidal cloud, then this approximation is likely to be in error by no more than 30%. 
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al 1995) and dense cores (Gaida, Ungerechts, & Winnewisser 

1984) mapped in molecular lines, as well as PPCs mapped in 
sub-mm dust continuum (Casali 1986; Ward- Thompson et al 
1994; Motte et al 1998). In contrast, radial profiles of cores 
with embedded stars lack such an extended, flat inner region, 
and tend to be more compact than PPCs (Motte et al 1998). 

An additional, theoretical, argument against the SIS model 
is that there is no known way for a cloud to evolve to a static 
SIS, since isothermal spheres are unstable to gravitational col- 
lapse for center-to-surface density ratios greater than 14 (Ebert 
1955; Bonnor 1956). Models of the collapse of isothermal 
clouds (Hunter 1977; Foster & Chevalier 1993) and of magne- 
tized clouds (Safier, McKee, & Stabler 1997) show fliat a l/r^ 
structure does develop during the coUapse, but it is not static; as 
a result, the accretion rate at early times is greater than that pre- 
dicted by the SIS model. However, at late times, both solutions 
behave similarly to the SIS. 

These considerations argue both against a singular density 
profile toward the centers of PPCs, and against the existence 
of a purely thermal region filling most of the PPC. On these 
grounds, we conclude that PPCs as a whole are unlikely to 
bear any close relation to (singular or non-singular) isothermal 
spheres. 

1.4. Previous Models of Self -Gravitating Cloud Equilibria 

Any complete model of a molecular cloud must include a de- 
tailed treatment of thermal balance throughout the gas, which 
entails a full description of the physical and chemical proper- 
ties at every point. While progress has been made along these 
lines (de Jong, Dalgamo, & Boland 1980; Falgarone & Puget 

1985) , polytropic models, which postulate a power-law relation 
between pressure and density at each point, 

P = K„p^^=Kpp'^K (1) 

are often adopted as a simpUfying assumption. Here, Kp is con- 
stant and jp and n are different versions of the polytropic index. 
Yet, spherical, non-singular, polytropic models of dense cores 
and Bok globules have met with only hmited success. Some- 
what surprisingly, the dominant problem lies not with details 
of chemistry, but rather in simultaneously reproducing the bulk 
properties — mass, radius, and density contrast — of the objects. 

Positive-index polytropes (n > 0, 7p > 1) cannot account for 
the observed variation of fine width with radius in molecular 
clouds since their temperatures decrease with radius (Shu et al 
1972). Gas temperatures in cores and globules are found to be 
either uniform, at T ~ 10 K (perhaps a factor of 2-3 higher 
in more massive objects and in regions with recent nearby star 
formation), or outwardly increasing (Bachiller et al 1987). Fur- 
thermore, on nearly all spatial scales, the nonthermal linewidth 
increases with size (Larson 1981). Negative-index polytropes 
(n<-l,0<7p<l; hereafter, NIPs), which have a tempera- 
ture T ocP/poc p^l" that increases with linear scale, can quali- 
tatively reproduce this behavior if the polytropic temperature is 
taken to represent the total velocity dispersion (Maloney 1988); 
i.e., Tffip(r) oc P{r)/p(r) cx o'^(r). The limiting case of the NIP 
as 7p ^ is the "logatrope", introduced by Lizano & Shu 
(1989) as a phenomenological model for turbulent pressure and 
studied further by McLaughhn & Pudritz (1996) and Gehman 
etal(1996). 

Molecular clouds and Bok globules are supported by sev- 
eral different pressure components — in addition to the thermal 
and turbulent pressures discussed above, there is also the static 



magnetic field. Mouschovias (1976) and Tomisaka, Ikeuchi, 
& Nakamura (1988) developed models of axisymmetric clouds 
supported by static magnetic fields and isothermal gas. Lizano 
and Shu (1989) extended these calculations to include turbu- 
lent pressure, modeled as a logatrope. More recently, McKee & 
Holliman (1999; hereafter MH) developed multi-pressure poly- 
trope models for clouds, in which each of the important pres- 
sure components is represented by an appropriate polytrope. 
Comparison of these models with observations was carried out 
by Holliman (1995). A phenomenological alternative to multi- 
pressure polytropes is the "TNT" model developed by Myers & 
Fuller (1992) and Caselh & Myers (1995). In this model, the 
density is assumed to be the sum of two terms, 

where the first represents the effect of thermal pressure and the 
second the effect of nonthermal pressure. By inserting this ex- 
pression into the equation of hydrostatic equilibrium, it is pos- 
sible to infer both the density at ro and the variation of the line 
width with r; reasonably good fits are obtained for most cores 
by an appropriate choice of the parameters ro and p. 

Fits of the envelopes of globules in the Clemens & Barvai- 
nis sample to a volume density law of the form n(r) oc r"*f 
indicate a rather narrow range: I <^ kp <^2 (Yun & Clemens 
1991). This further supports the use of NIPs, whose behav- 
ior asymptotically approaches a power law with precisely the 
same range of indices. [The singular polytropic sphere has 
n(r) cx r"^/P-7p) = ^-2n/(«-i) j Similar fits to dense core data pro- 
duce the same range of kp (Cernicharo, Bachiller, & Duvert 
1985; Bachiller et al 1987; Stuwe 1990). The use of NIPs in 
the thermal modeling of molecular clouds has been validated 
by de Jong et al (1980), Boland & de Jong (1984), and Falgar- 
one & Puget (1985), each of whom attempted to determine the 
equation of state relevant to dense molecular gas using detailed 
thermal and chemical balance calculations. While the repre- 
sentation of turbulence as a spatially varying isotropic pressure 
is clearly simplistic, Alfven waves, a significant contributor to 
MHD turbulence, have a polytropic index -jp = 1/2, correspond- 
ing to n = -2 (Walen 1944; McKee & Zweibel 1995). 

Non-singular NIPs that are truncated by a medium of con- 
stant pressure have well-defined stability properties. In keep- 
ing with the above discussion of virial equilibrium, it is of in- 
terest to compare the theoretical properties of critically stable 
NIPs with the observed masses, radii, and density contrasts of 
PPCs. Under the usual assumptions that the NIP is isentropic 
(see §2. 1), the mass, radius, and density contrast cannot exceed 
tiie values (Viala & Horedt 1974), 

where Mbe, Rbe, and {pd Ps)be are the mass, radius and den- 
sity contrast of the critically stable, bounded isothermal sphere, 
also known as the Bonnor- Ebert sphere (Ebert 1955; Bonnor 
1956). That is, the Bonnor-Ebert sphere (which has n = ±oo, 
7p = 1) sets the upper Umits of the physical parameters for 
any NIP, provided that the sound speed, c = {P/ p)^!^, in the 
isothermal sphere is the same as the surface value, c,, in the 
NIP. The same limits apply if Cg is replaced by the rms sound 
speed in the cloud, {c^Y'^, which is often more easily observed 
than Cs (MH). The fact that a 1 for observed PPCs ensures 
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that their masses and radii are in reasonable accord with the 
above, provided that c includes both thermal and nonthermal 
motions. However, the inferred density contrast between the 
intensity peak of many cores and the boundaries of their sur- 
rounding clumps (as determined from '^CO) is typically > 30 
(Stutzki & Guesten 1990; WilHams, Blitz, & Stark 1995rAn- 
dre et al 1996), which significantly exceeds {pd Ps)be = 14.04. 
The problem may be even more pronounced for Bok globules, 
which are expected to have a lower boundary pressure than 
cores since they are not embedded in GMCs. As stated suc- 
cinctly by Hasegawa (1988), who attempted to model the Bok 
globules L134 and L183, "...a single hydrostatic gas sphere 
cannot explain the high density contrast between the environ- 
ment and the center of globules." 

Attempts have been made to elude this difficulty. In their 
model of a quiescent Bok globule, Dickman & Clemens (1983) 
proposed a new stability criterion for NIPs; namely, that the 
central temperature, not the entropy, remains constant during a 
perturbation (see Appendix). This criterion was later adopted 
by Maloney (1988), who explicitly showed that it renders all 
bounded NIPs unconditionally stable. Of course, this removes 
the upper bound on the density contrast given above, allowing 
it to be fit to observations. Predictably, however, this can result 
in an overestimate of the mass, radius, and visual extinction 
of the object (Hasegawa 1988; Turner et al 1992). An impor- 
tant exacerbating factor in the modelling difficulties faced by 
both Dickman & Clemens and Hasegawa was their assumption 
that the total boundary pressure was equal to the thermal pres- 
sure (fth 3700 K cm"^) only. In contrast, current estimates 
including the interstellar magnetic field (but not cosmic rays, 
which pervade both the cloud and the ambient medium) give 
F ~ 2 X 10"* K cm"^^ (Boulares & Cox 1990). While this lessens 
somewhat the discrepancies they uncovered, it does not solve 
the problem. Nor does it change the fact that NIPs of any ra- 
dius are unable to fit the detailed shape of the line width-size 
relation (Paper II). 

Implicit in the stability analyses described so far is the as- 
sumption that the gas sphere responds isentropically, i.e. with 
spatially and temporally constant entropy, to perturbations ap- 
plied at its surface. While this is appropriate for an isother- 
mal gas alone, observations suggest a more complex situation, 
in which several distinct pressure components — primarily ther- 
mal, magnetic, and turbulent — coexist in the gas (Myers & 
Goodman 1988). As shown recently by MH, this necessitates 
a non-isentwpic description of cloud structure (the significance 
of this will be discussed in §2 below). Non-isentropic poly- 
tropes that accurately reflect the relative pressures of these dis- 
tinct components tend to be more stable than their isentropic 
counterparts. This improved treatment obviates the perceived 
need for the ad hoc, constant central-temperature restriction in- 
troduced by the above authors (see Appendix). 

1 .5. Motivation and Evidence for Core-Envelope Structure 

In this paper, we construct a model for PPCs that explicitly 
reflects their observed velocity structure; i.e., thermal and sub- 
sonic nonthermal motions on small scales (^ 0.1 pc), and su- 
personic nonthermal motions on larger scales. Theoretically, 
the transition between the two regimes is believed to occur via 
the damping of Alfven waves in regions of low ionization; i.e. 
near the density maximum where the optical depth to incoming 
FUV photons is largest. The critical wavelength for damping is 
of order 0.1 pc, under typical PPC conditions (§2.2). This sug- 
gests, in the apt description of Goodman et al. (1998), a view of 



PPCs as "islands of calm in a turbulent sea." 

Another theoretical argument for such a division comes from 
the theory of magnetized molecular clouds in the absence of 
wave motions. Central to the quasi-static picture is the notion 
of a magnetically-supported cloud, slowly contracting via am- 
bipolar diffusion (Mouschovias 1991). Over time, magnetic 
support in the central region decreases as the mass and den- 
sity increase there, while in the envelope, little changes. This 
leads naturally to a division between core and envelope, in that 
the latter remains "magnetically subcritical," while the core be- 
comes "magnetically supercritical" (McKee et al 1993). Even- 
tually, the core alone collapses, leaving behind the magneti- 
cally subcritical envelope that accretes more slowly (Safier et al 
1997). Inclusion of the effects of hydromagnetic waves should 
moderate the inward flow of mass via ambipolar diffusion, but 
the decay of the waves in the central region should allow it to 
operate as predicted, preserving — and perhaps enhancing — the 
distinction between core and envelope (see also Mouschovias 
1991). 

Studies of the thermal structure of clouds exposed to the in- 
terstellar radiation field also suggest a core-envelope structure. 
Falgarone & Puget (1985) studied the effect of external infrared 
and UV radiation on thermally-supported self-gravitating con- 
densations of gas and dust. For the conditions they assumed, 
they found equilibria featuring isothermal, molecular cores (of 
typical radius 0.1 pc and mass 2 Mq) surrounded by warmer 
and more diffuse, partially atomic, envelopes (with r ^ 1 pc, 
M ~ 10 Mq, depending on the UV field intensity and external 
pressure). More recently, detailed numerical calculations by 
Nelson & Langer (1997) demonstrated that the thermal evolu- 
tion of a 100-400Mq, uniform density cloud subjected to the 
interstellar radiation field leads to a core-envelope structure, in 
which the core region has roughly the same properties as those 
observed for PPCs in C'**0, while containing only one-fifth of 
the total mass. Elmegreen (1989) calculated the relative mass 
and radius fractions of molecular cores to atomic envelopes for 
spherical polytropes exposed to the interstellar radiation field. 
While this molecular/atomic core-envelope structure is gener- 
ally not relevant to PPCs, which are shielded by the ambient 
molecular gas, it may be relevant to the more massive Bok glob- 
ules, which have bulk properties quite similar to those predicted 
by the models of Nelson & Langer (1997), except for the large 
nonthermal linewidths in their envelopes. 

On the observational side, Bachiller & Cernicharo (1984) 
noted marked differences in the density and velocity structure 
of the large dark cloud B 1 depending on the molecular tracer 
used, prompting them to divide the cloud into an NH3 core 
(r ^ 0.3 pc, M < 60 M0), a C"*0 and '^CO envelope (r < 1 pc, 
M ^ 350 Mq), and a '^CO halo (r - 2 pc, M - 400 Mq). In- 
terestingly, despite these differences, the three regions showed 
negligible variations in kinetic temperature {T = 12±2 K). 
While this layering could simply reflect the density threshold 
at which different molecular species are excited, a later study 
of similar clouds in the Taurus-Auriga-Perseus complex using 
dust extinctions to derive H2 column densities found similar re- 
sults (Cernicharo et al 1985). On smaller scales. Fuller (1989) 
proposed the same division on the basis of higher-resolution 
NH3 and C'^O observations. 

As a particular example of a cloud one might keep in mind 
for our models, consider R65, a well-studied, starless clump in 
the p Oph complex (Loren 1989). '^CO maps of this clump give 
the foUowing estimates: T w 20 K, FWHM size « 0.77 x 0.44 
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pc, hh^ w 3.5 X 10^ cm-^ a w 0.39 km s"', and M w 21 Mq; 
no IRAS sources were detected within the FWHM contour. 
Note that we have assumed a H2 to '^CO abundance ratio of 
4.8 X 10^, in accord with the downward revision in Loren's 
original estimate by Bertoldi & McKee (1992). Using these 
values, and taking for R one-half of the geometric mean of the 
cloud dimensions, one finds a w 2.4. The dense core of R65, 
L1689B, has been detected in C'^O, NH3 (Benson & Myers 
1989) and, more recently, 1.3 mm dust continuum (Andre et al 
1996). The latter study gives: T w 18 K, FWHM geometric 
mean radius w 0.046 pc, w4.1 x 10'* cm"-', andMw 1 Mq. 
According to Andre et al, the derived density contrast between 
the intensity peak of L1689B and the edge of R65 is ^ 30, a 
figure much higher than simple equilibrium NIP models can 
accomodate, even at their critical states (i.e. pd ps < 14.04). 

In §2, we develop the theory of polytropes with core- 
envelope structures in detail, including an analysis of their sta- 
bility. After reviewing the properties of single-component poly- 
tropes in §3, we examine the physical properties of composite 
polytropes in §4. Some of the broader implications of these 
properties for star formation in PPCs are mentioned in §5. In a 
subsequent paper, we compare the models to available observa- 
tions. 

2. COMPOSITE POLYTROPES: GENERAL THEORY 

2.1. Locally Adiabatic Pressure Components 

Polytropic models have long been used to model molecular 
clouds. The discussion of the data on molecular cloud cores 
and Bok globules above suggests that two generalizations of 
the standard polytropic theory are needed in order to success- 
fully model these objects: First, as mentioned above, it is nec- 
essary to include several different pressure components, each 
with pressure Pj(r), so that the total pressure is P = Y^Pj- If 
7p y is the polytropic index for the jth pressure component, then 
the polytropic index for the resulting multi-pressure polytrope 
is given by (MH) 

where is the number of distinct pressure components. To 
simplify the analysis, we make the approximation in this paper 
that the spatially varying polytropic index "fpir) can be replaced 
by a suitable average, 7^ . 

Second, it is necessary to consider two (or more) spa- 
tially distinct polytropic regions, each with its own polytropic 
index — a composite polytrope. For a core and an envelope, we 
write 

^ core ^^-coreA^ ^ env '^^envA^ i 

where Kcose and K^w are constants. Composite polytropes have 
found wide application in stellar structure: Milne (1930, 1932) 
constructed such models for stars with degenerate cores. Hen- 
rich & Chandrasekhar (1941) considered models of a star with 
an isothermal core and an « = 3 (radiative) envelope, while 
Schonberg & Chandrasekhar (1942) investigated the effect of 
compositional gradients in such models. The radial stability of 
a composite polytrope with an isothermal core and an n = 3/2 
(7p = 5/3) envelope was studied by Yabushita (1975) in the con- 
text of neutron stars. More recently. Beech (1988) constructed 

'MH also consider the case in which the timescale for internal heat transfer is 
adiabatic." 



an analytic model for low-mass stars, joining an approximate 
solution for a n = 3 (7^ = 4/3) polytrope onto a n = 1 (7^ = 2) 
envelope. 

The specification of 7^ in each of the two regions suffices for 
the determination of an entire equilibrium sequence; e.g., the 
total mass as a function of central density. In order to decide 
the issue of stability for members of that sequence, however, 
one needs to know the response of each pressure component 
to small perturbations. We shall assume in this work that: (1) 
the perturbations are adiabatic, and; (2) for each component, 
the time scale for internal heat transfer — a direct consequence 
of the perturbation — is long compared to the dynamical time 
scale. Following MH, we refer to the latter as "locally adi- 
abatic" pressure components^. Such components satisfy the 
usual adiabatic relation (Ledoux & Walraven 1958), 

5 In Pj = -fj6 In p, 

where the different pressure components are assumed to be 
thermally insulated from each other. The adiabatic indices jj, 
then, describe temporal variations of the pressure components, 
while the polytropic indices, ^pj, describe spatial variations. 
For locally adiabatic polytropes, MH show that the overall adi- 
abatic index is 




For the purposes of this discussion, we shall assume that, as in 
the case of the polytropic indexes, this spatially varying adia- 
batic index can be approximated by a mean adiabatic index in 
each region, 7. Calculations using multi-pressure polytropes 
(MH) do not suggest a simple prescription for determining 7^ 
or 7 from the properties of the constituent polytropes. However, 
these calculations do indicate that 7^ and 7 are weighted means 
in the sense that if there are two pressure components with adi- 
abatic indexes 71 and 72, for example, then 71 < 7 < 72. In 
Paper 11 we shall infer values of 7^ and 7 from observations. 

When 7p = 7 in a region, the gas is isentropic there (i.e., it 
has a constant specific entropy), and as a result no internal heat 
transfer is required to maintain the polytropic structure after a 
perturbation. Most studies of polytropic models for interstel- 
lar clouds have assumed that the clouds are isentropic. More 
generally, however, 7^ y^j. For example, small-amplitude, un- 
damped Alfven waves have ^p,„ = 1/2 and 7,, =3/2 (McKee 
& Zweibel 1995). Note that when 7 < 7^, polytropic spheres 
are convectively unstable according to the Schwarzschild crite- 
rion (Cowling 1941); none of the known pressure components 
are in this regime, however. For the pressure components in 
molecular clouds, we have 7th = Jp,th for the isothermal gas; 
7b > 7p,B for the magnetic field in order that the cloud be stable 
against the interchange instability; and, as we have just seen, 
7w > 7/),w for the Alfven waves. Since each pressure compo- 
nent has 7 > jp, we expect 7 > 7^ in each part of a composite 
polytrope model for molecular clouds. This is significant, since 
non-isentropic polytropes are more stable than their isentropic 
counterparts (MH). This is essentially due to the fact that as 7 
approaches 4/3, single-component polytropes are uncondition- 
ally stable (Ledoux & Walraven 1958). 

short compared to the dynamical timescale; they term such components "globally 
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2.2. Core Properties 

All available evidence suggests that the dense, thermally- 
dominated core of a PPC is an unusually quiescent region, in 
comparison to its surroundings. We presume, as have others, 
that there is a scale above which the cloud is supported against 
its own weight primarily by a mean magnetic field and MHD 
waves (Shu, Adams, & Lizano 1987; McKee et al 1993); for 
typical low-mass star-forming regions in Galactic molecular 
clouds, this length scale is about 0.1 pc. In the core, how- 
ever, the observed near-thermal velocity linewidths suggest that 
the waves important for maintaining equilibrium on large scales 
are damped at smaller scales; specific damping mechanisms for 
MHD waves were discussed by Zweibel & Josafatsson (1983). 
The cutoff wavelength for damping along the direction of the 
field has been estimated as (Kulsrud & Pearce 1969) 



Acut = 0.13pc 



B 



V lO^^cm- 



Cc 



-1/2 



(4) 



where B and n^^ are characteristic values of the magnetic field 
strength and H2 density in the relevant region, and (^cr is the 
cosmic-ray ionization rate. We have assumed the standard bal- 
ance between ionization and recombination in order to compute 
the coefficient (McKee et al 1993). However, (cr can range be- 
tween 10"'^ CcR < 10"'^ depending on the degree of shield- 
ing (Caselli et al 1998), meaning that Acm is uncertain by as 
much as a factor of 3 in either direction. Nevertheless, equa- 
tion (0) suggests an important constraint on the models; i.e. 
that the thermally dominated core region should have a radius 

^core — Acut/2. 

We now turn to the poly tropic and adiabatic indices. By con- 
struction, jpxh = 7th = 1 ■ The value of 7^.^ is somewhat more 
uncertain; it depends both on the degree of isotropy in the con- 
traction of the cloud from a more uniform initial state, and on 
the amount of ambipolar diffusion that has taken place. For 
a PPC formed by compressing a uniform density cloud under 
strict flux conservation (i.e. no ambipolar diffusion), 7^,.^ = 4/3, 
while for a PPC of "infinite age" (i.e. maximal ambipolar dif- 
fusion), jp,B = 0. As an intermediate value, we shall assume 
Jp.B = 1 in the core of our composite polytrope, indicating a 
moderate amount of ambipolar diffusion. Since this is identi- 
cal to the value of 7^ for the gas pressure, the other pressure 
component in the core, we adopt 7^ core = 1 ^ thus, the isother- 
mal solution can be used for the core. In a subsequent paper, 
we shall use observations to estimate a range of likely values 
for 7pcore' show that this too is consistent with the use 
of the isothermal approximation. Since the redistribution of 
flux is negligible on the time scale of the perturbation, we take 
7b = 4/3, appropriate for adiabatic, flux-conserving compres- 
sion. We expect the adiabatic index for the core to have a value 
intermediate between that of the thermal gas and the magnetic 
field: l^S7core;$4/3. 



and of wavelengths shorter than variations in the background 
density, they exert an isotropic pressure obeying P cx p'*'' with 
7p =1/2 (this relation was originally derived by Walen 1944). 
Further, the adiabatic index of such waves is 7,,, = 3/2. Now, 
were Alfven waves treated as a locally adiabatic pressure com- 
ponent, they would be unconditionally stable, since 7„, > 4/3. 
However, as argued by MH, undamped Alfven waves are more 
likely to behave as a globally adiabatic pressure component, 
since wave energy tends to flow from high to low density re- 
gions in a centrally-concentrated cloud. These authors demon- 
strated that single-component, globally adiabatic spheres are 
generally less stable than their locally adiabatic cousins. In par- 
ticular, a critically stable cloud with 7,, = 1/2 and 7,,, = 3/2 has 
a finite center-to-surface pressure contrast of 4. 15, and a density 
contrast of 17.26. 

In our description, in which Alfven waves are treated in an 
averaged sense, we make the approximation that the waves, 
like the mean magnetic field, are locally adiabatic. As shown 
by MH, the location of the critical point for a multi-pressure 
polytrope including Alfven waves is the same as that for a 
similar multi-pressure polytrope in which the Alfven waves 
are replaced by a locally adiabatic pressure component with 
"fpj =1/2 and with a suitable value of -fj. For a cloud supported 
entirely by Alfven waves, the value of 7 for the equivalent lo- 
cally adiabatic component is 7„,eq = 0.6056. However, as other 
pressure components become important, the value of 7„._eq in- 
creases, which complicates the choice a suitable value for 7. 
In any case, since 7 is some weighted mean of the 7/s of the 
individual components, we expect 7 < 4/3. 

We now proceed to the actual construction of composite 
polytropes, starting with the illustrative example of the com- 
posite isothermal sphere. Note that in what follows, we omit 
the overline indicating the averaged character of 7^ and 7; for 
the theory, it matters only that they are constants. 

2.4. Isothermal Cores and Envelopes 

In order to gain a feel for composite equilibria, we begin with 
the illustrative example of the composite isothermal sphere. 
Consider a core of isothermal gas at temperature T^oTe sur- 
rounded by an isothermal envelope at temperature Tenv > Tcoie- 
Let Ccore, env denote the constant isothermal sound speed in the 
core and envelope, respectively. Define the jump in conditions 
at the interface by 



^env \ ^nv -, 
Ti=\ I = — > 1, 



env 

*^core / ^ core 



(5) 



where in the second step we have assumed both regions to be 
of the same composition. 

Introducing the customary dimensionless variables ^ and 
V'(0 (Chandrasekhar 1939), the density, radius, and mass 
within the core are given by 



2.3. Envelope Properties 

The observed increase of nonthermal linewidth with size to 
supersonic values in molecular clouds strongly suggests the 
presence of MHD turbulence at scales larger than about 0. 1 pc. 
We therefore include a third, Alfven-wave pressure component, 
Pk, in the total mean pressure of the envelope. 

The assumption of a polytropic relation for Alfven waves 
was examined by McKee & Zweibel (1995). They showed 
that as long as the waves are hnear in amplitude, undamped. 



p= PcSXpi-lp), 



{AnGpcfl^ 
Mir) = 47rp,rie^, 



(6) 
(7) 

(8) 



where pc is the density at the center (r = 0) and r^^ is a charac- 
teristic radius proportional to the Jeans length in the core. In the 
interest of notational simplicity, it is desirable to retain ^ and ip 
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in the envelope as well as in the core. In the envelope then, 



p = p^exp(-V'), 

n 



(9) 



(4^Gp,)i/2 



M{r) = 47: p.r^e^, 



where pe is an undetermined constant. 

At the interface of core and envelope, the radius, pressure, 
and mass must be continuous (but not the density; the present 
class of models differs from composite stellar models in pre- 
cisely this sense). The above expressions therefore give 



CcorePceXp(-(/'-) = C^„^p,eXp(-0+), 



(10) 



'-envPe r+J 



where a prime denotes the derivative of a function with respect 
to its lone argument, a subscript '-' indicates a quantity evalu- 
ated on the core side of the interface, and a subscript '+' indi- 
cates a quantity evaluated on the envelope side of the interface. 

The equilibrium in both regions is governed by the Lane- 
Emden equation. 



^(eV)' = exp(-V'). 



(11) 



The unique solution of equation (|1 1[), subject to the standard 
(i.e. regular) boundary conditions, 

^ = 0, V' = at C = 0, 

is the isothermal function tabulated by Chandrasekhar & Wares 
(1949). Thus, for a given dimensionless core radius, = 
fi/rch, ip-, and "0- known (a subscript / denotes quantities 
that are continuous across the interface). Equations ( p^ are 
then three equations in three unknowns: i/;^, and Pe/pc- 
Solving these, one obtains the following relations: 

1p+ = tp- + 2lnTi, lll'+=ll^'_/Ti, Pelpc = Ti. (12) 

Equations ( p^ give starting values for the integration of 



equation ( |1 1| ) out to some radius r, = rchCs^ where an outer 
boundary condition is imposed. We imagine that the compos- 
ite cloud is ultimately surrounded by a tenuous (compared to 
the envelope) intercloud medium, whose entire dynamical ef- 
fect is to exert a finite, constant pressure everywhere on the 
outer boundary. This situation is effectively the same as is as- 
sumed for the well-known bounded isothermal sphere (Ebert 
1955; Bonnor 1956) except that we suppose Ps to be comprised 
of both thermal and nonthermal components. Using equations 
(||), (||), (^) and (|l2|), we find that this pressure is related to the 
central pressure, Pc = c^okPc, via 



■ = T; exp(^s). 



(13) 



Note that the corresponding density contrast is simply Ti{Pc/Ps)- 
Alternatively, one may specify a given configuration by its ratio 
of core-to-total radius and/or mass. The definitions of M and r 
along with equations ( p^ yield 

Mi _ __i ($V0- 



(14) 



Note that when expressing products of quantities, any one of 
which is not continuous across the interface, we retain the '-' 
subscripts for the entire product. 

2.5. Isothermal Cores with Polytropic Envelopes 

We next consider composite poly tropes with isothermal cores 
and polytropic envelopes as an approximate model for molec- 
ular clouds. (In principle, the core of a composite polytrope 
can have an arbitrary polytropic index, but we restrict consid- 
eration here to the special case in which it is isothermal.) As 
remarked above, composite polytropes have long been used in 
studies of stellar structure. The general matching conditions for 
two-component spheres of differing polytropic indices n were 
presented by Chandrasekhar (1939). 

The Lane-Emden equation governing the envelope is 



i»+ii 

n+ 1 



= ±1, 



(15) 



where the suitably non-dimensionalized variables for the enve- 
lope are 

p= PeO , P = KenwP ' = AenvPg ' & , 



l»+l| 

4ttG 



-K„ 



1/n-l 



1/2 



M{r) = At: /\p,rl^eO' 



where we have used the continuity of ^ to express the mass in 
terms of characteristic radius for the isothermal core, rch. Re- 
call that n and 7^ are related via 7^ = 1 + 1 /n. The matching 
conditions at the core-envelope interface are 



CcorePceXp(-(/'-) = K^nyP, 



-1/2. 



^corePc 

Pcipi = Ap,e'^ 



1 1 I l/«-Kl/2 



(16) 



Generalizing the definition of r,- (eq. ^ to the case of a poly- 
tropic envelope, we have 



' ~ P+ Pe01 

The first two of equations (|6|) can be combined to give 



(17) 



exp(ij-)ei^^ = \n+l\. 



(18) 



Equations (p^-(p^ then yield the jump conditions: 

0+ = [Tf |n+l|exp(^_)]'/<'-«', 

Pel Pc = T~^e.xp{-%l;-)ei". 



(19) 



These are the equivalents of equations (jl^) in the case of a poly- 
tropic envelope. 
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The ratio of core-to-total radius is again given by equation 
(p4|), while the mass ratio is 

Finally, the center-to-surface pressure and density contrasts are 
given by 



^ = |n+l|[r, exp(V'-) W 



\n-\2f\-n—\ 



(20) 



As discussed in §§2.2-2.3 above, the adiabatic indexes for 
both core and envelope are expected to be less than 4/3. As a 
result, the composite polytrope must be bounded by a medium 
of finite and constant pressure at its surface (MH). Also note 
that the envelope temperature T ex P/ p cx p'^" is an increasing 
function of ^ in NIPs (Shu et al 1972; Viala & Horedt 1974). 
As mentioned in §^^, we interpret the departure from isother- 
mality in the envelope as arising entirely from the presence of 
mean magnetic fields and Alfvenic turbulence, approximated 
as isotropic pressure components. The latter appear as nonther- 
mal contributions to the total sound speed, given by (McKee & 
Zweibel 1995) 



C ('•) = --T- = CTth 

p{r) 



+ o^.v('")+o^l('"): 



(21) 



where crth is the (constant) thermal speed, cr„(r) is the ve- 
locity dispersion associated with Alfven waves, and VA{r) = 
B{r)/\/4-TTp{r) is the Alfven speed. In our model, (!„ is the 
same as the nonthermal velocity dispersion (T„t. Note that the 
total velocity dispersion is related to c(r) by 



a(r) = [ai + al(r)] i/2 = [ ^c2(r) + 1^ - 



3- - 3-th-3''A('-)] 



1/2 



(22) 



2.6. Stability of Composite Configurations 



While polytropic models of molecular clouds clearly repre- 
sent a simplification of the real situation, they have a distinct ad- 
vantage over more detailed models, in that the analysis of their 
stability is straightforward. We note, however, that previous 
studies disagree as to whether such a cloud responds isother- 
mally or adiabatically to radial perturbations applied at its sur- 
face. We review this situation in the Appendix, and justify our 
choice of the adiabatic stability criterion adopted here. 

The equation of radial motion for small, adiabatic perturba- 
tions about equilibrium is, quite generally (Eddington 1926; 
Ledoux & Walraven 1958) 



4 1 t/7 GMir)p 
r J dr r^P 



uP'p _^ 3 d^ 



jP jr dr 



dh 
dr 

4\ GM{r)p 

1 



r^P 



h = 0, (23) 



where h = 6r/r is the relative displacement of a fluid element at 
a radius r and ut is the frequency of the oscillations. We assume 
that the adiabatic index, 7, can take on any (positive and con- 
stant) value in either the core or the envelope. Thus, in general, 
7 has a discontinuity at the interface. For the reasons outlined 



§ 1.2, we confine our interest to critically stable clouds, setting 



: in what follows. 



Under these assumptions, equation 
isothermal core. 



becomes, in the 



^ ^core 4 

~d^ 



dhr. 



d^ 



-(3-4/7coie)-^/?core = (isothcrmal). 

(24) 

When 7coi-e = 1, Yabushita (1974) showed that there exists an 
exact solution of equation (|24|); namely. 



exp(i/)) - 1 



where e is a small, dimensionless constant. 

The solution of equation ( p4| ) must be matched onto the cor- 
responding solution in the envelope, h = h^n- For an isothermal 
envelope, the governing equation is identical in form to equa- 
tion (p4[); for a polytropic envelope, one has 



4 e' 



dhf. 



di 



+ (3 -4/7e„v)(n+ 1)— /Jenv = (polytropic). (25) 

By definition, h must be continuous at the interface; i.e. 
hemi^i) = hcore(£,i)- Its derivative must satisfy the condition 
(Ledoux & Walraven 1958) 



dhf. 



di 



Tcore / dh^, 



7env V "^f 



3/Zc 



6' 



7c. 



7ei 



which takes into account the possibility of a discontinuity in 7 
across the interface. 

The condition for the critical point is that the Lagrangian 
variation in the pressure at the boundary vanish, {5P/5r)M, = 
(Ledoux & Walraven 1958; MH). This places the following re- 
striction on /ipn„: 



=0. 



(26) 



This condition holds for both isothermal and polytropic en- 
velopes, and differs from those expressed by equations (Q) and 
(^0|), in that it holds only for critically stable states. Finally, note 
that although we shall take « < -1 (0 < 7y, < 1) in what follows, 
the results of this section hold for all n ^-1. 

2.7. Solution Procedure 

In practice, the solution procedure for critical states is as fol- 
lows: 

(1) Choose a model by specifying five parameters: 

n (or 7p), 7core, 7env, and r,-. 

(2) Solve the Lane-Emden equation for an isother- 
mal sphere on the interval [0,^,]; i.e., find ?/;(4) and 

(3) Using V-, /!core(6). and (dh^om/dO- as 
starting values, calculate 9'^, henvi^d, and 
{dh^nw I d^)+. These are the starting values for the 
integration of the two ODES (ITTb and ( ^4| ) (isother- 
mal envelope) or (|5]) and ([25|) (polytropic enve- 
lope). 

(4) Integrate the system outward in ^ until the outer 
boundary condition (Eq) is satisfied. 
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The resulting unique solution yields, as dimensionless output 
parameters: (i) the center-to-surface pressure and density con- 
trasts; (u) the core-to-total radius ratio; and (iii) the core-to- 
total mass ratio. 

3. SINGLE-COMPONENT POLYTROPES: A REVIEW 

Consider the variation of the mass M of a single-component 
NIP (-00 < « < -1; < 7p < 1) if the surface pressure and en- 
tropy are held fixed while the central density is increased. For 
fixed surface conditions, the behavior of M follows that of the 
the dimensionless mass ^ defined by MH, 



M 



M 



(27) 



The sequence of n versus pd Ps ranges from zero at pd Ps = 1 
up to some maximum /imax, beyond which it oscillates about the 
value corresponding to a singular polytropic sphere as Pc/Ps — * 
00 (Chandrasekhar 1939; MH), i.e. 



Moo ~ 



1/2 



(4-37,)'/^-^y^ 



(2-7p)2 



1/2 



-3|'/>+lp/2 
(n-l)2 ■ 



The sequence is qualitatively similar for all values of < 7^ < 
6/5 : the maximum possible mass of a poly trope in this range is 
determined by the conditions at its surface. (The same conclu- 
sion holds if Cs is replaced by the rms value (c2)i/2.) 

The largest stable density drop along this sequence occurs at 
the critical mass p^, which is less than or equal to /imax- In 
the isentropic case, j = ^p, one has pa = Mmax (Shu et al 1972). 
In this case, all higher density, lower mass models are gravita- 
tionally unstable. When 7 > 7^, p^r < Pnax, even though the 
corresponding {pc/ Ps)a is greater than in the isentropic case. 
This behavior is illustrated in Figure 1, which displays the par- 
tial equilibrium sequences of both an n = -3 (7^ = 2/3) poly- 
trope and the isothermal sphere. Recall that when 7 < 7^, both 
polytropic and isothermal spheres are convectively unstable, ac- 
cording to the Schwarzschild criterion (Cowling 1941). 

Sometimes we know the conditions at the center of the poly- 
trope, and can ask how much mass can be supported. In 
this case, the equilibrium mass can increase without limit as 
the surface pressure and density are lowered (Viala & Horedt 
1974); however, such equilibria are stable only if they are suf- 
ficiently non-isentropic. Noting that Cj = Pj/p^ oc p]^", we have 

M cx ppf "^^^"; as ps/ Pc^O for fixed pc, p — > Poo and M — > 00 
for n > 3 and n < -1 (0 < 7p < 4/3). 



4. THE COMPOSITE ISOTHERMAL SPHERE 

4.1. Pressure and Density Structure 

In Figures la and b, we plot, for various values of r,, the 
critical density contrast of the composite isothermal sphere 
(CIS), (pc/ Ps)ci.cr, as a function of rj/r, and Mj/Ms, respec- 
tively. We first present results for the isentropic case only, i.e. 
7core = 7env = lp= 1 ■ Notc that the corresponding pressure con- 
trasts are simply a factor of r,- lower than {pd Ps)ci,ci- The left- 
hand border of Fig. 2 corresponds to the case of an almost pure 
envelope with a core at the origin, whereas the right-hand bor- 
der corresponds to the case of an almost pure core with an en- 
velope at the surface. In each case, the limiting values of the 
curves at the borders are {pd Ps)ci,ct Tiipc/Ps)BE, indicating 
that (Pc/Ps)ci,cv (Pc/Ps)be, as expected. Note that for values 



of Ti ^ 5, configurations with the same ri/rs can have differ- 
ent (pc/ Ps)ci.cr- These equilibria are still unique, however, as 
evidenced from the absence of this feature in Fig. 2b, which 
employs the Lagrangian variable M,- /A/,. 

For each value of t,-, there is a maximum in {pd Ps)ci.a, 
which increases with increasing r, . As r, — * oo, the density con- 
trast across the interface becomes infinite, and thus an infinite 
center-to-edge density contrast can be supported. The maxi- 
mum pressure contrast, however, is finite, and can be estimated 
as follows. Since the core becomes less important with increas- 
ing Ti (it is already only 5% of the total radius and 1% of the to- 
tal mass at r, = 5), one can view this region in the large-r, Umit 
as a tiny, cold and dense "pea" in the midst of an essentially 
infinite, hot and tenuous background. These are precisely the 
conditions pertaining to the bounded isothermal sphere. Thus, 
on the scale at which pressure gradients in the core are signif- 
icant, ^ <^ ^i, the maximum pressure contrast across the core 
alone will be (Pc/Pda ~ (Pc/Ps)be- On the other hand, on the 
scale at which pressure gradients in the envelope are signifi- 
cant, ^ ^i, the region within the core is essentially structure- 
less; its only significance for the equiUbrium is to set the inter- 
face values of various quantities. Thus, the expected maximum 
pressure contrast in the envelope alone will again approach the 
Boimor-Ebert value. Together, these results imply a total pres- 
sure contrast in the large-r, Umit of 



C/,cr 



~ 197. 



BE 



The critical stability curves of pressure contrast versus r,/rj 
do indeed approach this asymptotic limit for r,- ^ 1 (Fig. 2c). 
Correspondingly, the total pressure contrast for a composite 
polytrope consisting of A' components, each separated from the 
other by a large value of r,, approaches (14.04)^^. 

4.2. Radius and Mass 

The dimensionless outer radius of the critical CIS can be 
much larger than that of the Bonnor-Ebert sphere; e.g. when 
Ti = 2 and rj/r, = 0.1, it extends to ^c/,cr = 15.7, versus £_be = 
6.45. Note that £,ci,ci is a function of the core-to-total radius 
ratio r,7 r^. As this ratio approaches zero (i.e., a pure envelope), 
one finds that ^c/,cr — > Tj ^be, or ^c/,cr = 12.9 in the case consid- 
ered here. Hence, the outer radius attains a maximum at some 
intermediate ri/r^. 

Because the temperatures of the densest observed regions 
of molecular clouds consistently fall in the range 10-20 K, it 
is convenient to normalize the mass with respect to this value 
rather than the surface value. On the other hand, the pressure at 
the surface of a cloud is dictated by the ambient conditions, and 
is often better determined than the pressure in the cloud inte- 
rior. We therefore introduce a hybrid dimensionless mass based 
on the temperature in the core and the pressure at the surface of 
the envelope: 

M(r) 



m{r) -- 



diG'p.y/'' 



(28) 



For fixed Ccore and Ps, m(r) is directly proportional the dimen- 
sional mass M(r). By contrast, the dimensionless mass p intro- 
duced above is normalized to the characteristic Jeans mass at r; 
the two are related by 



l^(r) ^ Ccore 

m(r) c^(r) 



P(r) 



1/2 



1 



P(r) 



1/2 



(29) 
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where 



For a CIS, r = 1 in the core and r = r, in the envelope. 

The total dimensionless mass of a composite isothermal 
sphere is then 

?MC/ = 'Mcore + »Jenv,I, (30) 



where 



pl/2^3/2 



-Me, 

core 

'^^env.l — A ^^env.I ~ 



4-TT- 



= 47r- 



-1/2 



(31) 



-1/2 



r,[(?V).-(?V)J. 



In the limit that the CIS is a pure core (r,/rj — > 1, ^_ — > ^j), 
nici = Wcore must approach mgi, the dimensionless mass of the 
bounded isothermal sphere. 



mB,= —^—M= Utt- 



-1/2 



(CV) 



(32) 



The last step, which follows from equations (^ and (^), shows 
that mbi indeed agrees with mcore from equation (pT) in this 
limit. In the opposite limit in which the CIS is a pure en- 
velope (xi/rs ^0, C+ ^ 0), equation ( p9[ ) implies mci 
(c/ccoiefniBi = rfmBi- The maximum mass of the CIS occurs 
for the pure envelope (see Figure 3), and since niBi < niBE, we 
find that the upper limit on the mass of a CIS is 



ma < T; mbe- 



(33) 



The above expressions are also valid for the critically stable 
CIS. The critically stable pure core or envelope is the Bonnor- 
Ebert sphere with ^j,be = niBE = 1.1822. Note that because the 
critical values of Pi/P.,, Wcore, and menv.i depend on the frac- 
tional core radius rj/r,, inciter a function of ri/rg. We plot 
inci.a/inBE, in the isentropic case, for various values of t,- in 
Figures 3a and b. The individual core and envelope components 
are also shown. Note that in general, as r, increases, the core 
becomes a lesser contributor to the overall mass of the cloud. 
As r, — > oo, we recover the pure envelope case. Finally, note 
that the above mass expressions are quite general; they apply to 
both isentropic (7 = 7^,) and non-isentropic (7 ^ 7p) spheres. 

4.3. The Non-isentropic CIS 

Consider now the much larger class of non-isentropic, criti- 
cal CISs; i.e. those with 7 1. We a re m ost interested in the 
range 1 < 7 < 4/3, as discussed in §2.3. Single-component 
isothermal spheres with 7 in this range were examined by 
Yabushita (1968), who showed that the dimensionless critical 
radius and pressure contrast of such spheres approached infinity 
as 7 ^ 32/25. The critical radius and pressure contrast remain 
infinite for 32/25 < 7 < 4/3 (MH). 

These features of single-component spheres can largely be 
carried over to the CIS. That is, CISs with 7coie > 1 and/or 
7env > 1 have pressure contrasts in excess of those found in the 
isentropic case, but smaller masses. The maximum pressure 
contrasts (i.e., as a function of r,/rs) for CISs with t, = 2 and 
Ti = 3 are plotted as contours in the (7core, 7env) plane in Figures 
4fl and b. Note that as Jcok and 7env approach the value 32/25, 
the various physical quantities diverge, and it proves difficult to 



integrate the equations all the way up to the singular point. This 
is a feature of Fig. 4 and several of our subsequent plots; how- 
ever, we need not be concerned that any essential qualitative 
behavior is being missed. 

To outline the important qualitative changes in the non- 
isentropic CIS, let us restrict consideration to r, = 1 and 1 < 
7 ^ 1.25. Then the critical mass extrema always obtain in 
the limit of a pure core or envelope, depending on the rela- 
tion between 7coi-e and 7env If 7core > 7env, then the minimum 
(maximum) mass is given by a pure core (envelope); and con- 
versely if 7(;ore < 7env Further, these extrema are given by the 
usual bounded isothermal sphere sequence beyond the mass 
peak (Fig. 1). When t, > 1 and/or Jcok or 7env <; 1.25, the 
mass extrema no longer necessarily occur for a pure core or 
envelope: the values of rz/rj corresponding to the mass peak 
depend on the exact values of r,, 7core, and 7env. As previously 
discussed, 7env > 1 'Wenv.i < i^BE, cvcn in the pure envelope 
limit. On the other hand, for r, > 1, this effect competes with 
a relatively larger mass contribution from the envelope, due to 
its increased temperature (§4.2). This combination sometimes 
leads to a maximum in nici.a- (< ^ibe) at intermediate r,/rj, an 
effect manifested by the slight asymmetry in the (Pc/A)c/ct 
contours between corresponding pairs of (7core,7env) in Fig- 4fl. 

4.4. Summary — CIS 

Whereas the bounded isothermal sphere (usually assumed 
to be isentropic, so that 7=1) has a maximum density and 
pressure contrast between the center and the surface of a fac- 
tor 14.04, composite isothermal spheres can have much larger 
drops: for two components with a temperature jump r,, the 
pressure drop can be up to a factor (14.04)^= 197 (in the limit of 
large t,), and the density drop a factor r, larger. Non-isentropic 
isothermal spheres, whether composite or not, can have very 
large density and pressure drops provided 7 is close to 4/3. 
We attribute the increased stability of composite polytropes and 
non-isentropic polytropes to the fact that they both have an en- 
tropy that increases outwards. The entropy is proportional to 
ln(P/p''), which increases outward as (7- 7p)ln(p(./p) for a 
non-isentropic poly trope and as 7lnT, across a jump in a com- 
posite polytrope. 

The core mass of a CIS can never exceed its Bonnor-Ebert 
mass (mcoie £ t^BE = 1-18), and the mass of the entire CIS can 
never exceed the Bonnor-Ebert mass of the envelope (pci < 
1.18). As shown in Fig. 3 and in equation (^3|), however, the di- 
mensionless mass of the CIS can exceed mBE whenever r, > 1, 
since mci.a- — > ^fmBE in the limit of a pure envelope. The di- 
mensionless mass of a non-isentropic CIS is less than that of an 
isentropic one (with the same r,). 

In order to compare the mass of a CIS with that of an actual 
dense core or globule, we need to re-express equation ( |30| ) in di- 
mensional form. Scaling to the representative values T^ok = 10 
K and f ,/fc = 2 x 10"* cm"^^ K (Boulares & Cox 1990; we have 
subtracted the cosmic ray pressure, since cosmic rays penetrate 
the molecular gas as well), one finds 



CI 



^core 



-nici = l.l Mq 



lOK 



'2x lO^cm-^KX'/^ 
X I I mc/, 



P./k 



(34) 



where a mean molecular weight of 2.33 niH has been assumed. 
The quantity nici, scaled by niBE, can be read off from Fig. 
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3. Assuming these representative values, the mass of the isen- 
tropic, critical CIS falls within the range 2.6 Mq < Mci.a < 
2.6 t?Mq (see eq. |3^. The total molecular masses of dense 
cores and globules derived from ammonia and C"^0 measure- 
ments range between a few Mq and a few tens of Mq (§1.1), 
so critical CISs with r, < 5 are in good agreement on this count. 

Since observations show that PPCs and Bok globules both 
have substantial nonthermal motions in their envelopes, CISs 
are primarily of academic interest. We turn now to the more re- 
aUstic case in which the envelope has a substantial nonthermal 
pressure. 

5. THE COMPOSITE POLYTROPIC SPHERE 

In models having a polytropic envelope, we assume that the 
thermal velocity in the envelope is the same as that in the core, 
fth, env = fth, core- The polytropic temperature, T{r) (x c^(r), 
may be interpreted as the sum of the kinetic temperature, Tkin 
(fixed at the core value, Tcoie), and an "effective" temperature, 
Tesir) oc (3o'^ + v^)/2 (see eq. ^1]). At the inner edge of the en- 
velope, reff.+ = (r, - 1) Tcoie, so that when t, > 1, the jump in T 
across the interface is assumed to be entirely due to a nonzero 
value of Teff(ri).^ We focus on the case of NIPs for the envelope 
(n<-l, 0<7;, <1). 

5.1. Pressure and Density Structure 

Beginning again with the isentropic case, in Figure 5fl we 
plot, for various negative values of n and for t, = 1, the criti- 
cal density and pressure contrasts of the composite polytropic 
sphere (CPS). In Figs. 5b and c, results are given for t, = 2 and 
r, = 3, respectively. For t, = 1, these curves approach the cor- 
rect values for pure polytropes at rj/r, = (Shu et al 1972)^ 
and the Bonnor-Ebert value at ri/r,, = 1. For t, > 1, the density 
curves approach t, times these same limiting values. Observe 
that CPSs attain peak critical densities that are typically a few 
times larger than their pure polytropic counterparts for t, = 1, 
and up to an order of magnitude (or more) larger for r, > 2. 
Also note that for r, > 1 and \n\ sufficiently larger than 1 (i.e., 
jp ^1), (Pc/Ps)a can exceed {Pc/Ps)be for a finite range of 
r,7 Tj, just as in the case of the CIS (see Fig. 2c). 

Fig. 5 also shows that the limit n ^ -I {jp 0) is quite 
anomalous, in that the density contrast (but not the pressure 
contrast) greatly exceeds that of other CPSs over a large range 
of r,/rj. This can be understood as follows: For small r,/rj, the 
CPS is mostly constant-pressure envelope, and is thus highly 
unstable. But as ri/r,, increases, stability is recovered due to the 
presence of the isothermal core. More importantly, the high 
compressibility of the n ~ -1 envelope has two distinct and 
complementary effects: first, the core is relatively insensitive 
to perturbations occurring at the cloud surface; and second, a 
relatively small pressure drop across the envelope corresponds 
to a large density drop. As in previous cases, these effects are 
only enhanced when r, > 1 . 

5.2. Radius and Mass 

As in the case of the CIS, the dimensionless outer radius of 
the critical CPS has a maximum at an intermediate value of 
ri/rs- This maximum is always greater than the Bonnor-Ebert 
critical radius, and therefore also greater than S^pct, the critical 
radius of the pure polytrope (Shu et al 1972). As r,/ rj ^ (i.e., 

*In some cases, e.g., Bok globules exposed to the interstellar radiation field, 
and TefKr) = T(r)-Ti,j^{r). 

'Note, however, the following error in Shu et al's Table 4. While (Pc/Ps)a — > 



a pure envelope), one finds that ^cp.cr ^ t-^""""''/^" 
e.g., ^cp,cr = 6.77 when n = -3 and r,- = 2. 

The mass of the CPS, in dimensionless form, is 

OTc/> = OTcoie + '«env,P, 



where mcore is given by equation (|3 

PI'^G^I^ ( Pr 

menM,P = 7 Menv.P = { 4-7r- 



-1/2 



At, 



\n+l\ 



X [(eo%-(eo%i 

and T, = T(rs) = c\l c\^^^. In the limit that the CPS is a pure core, 
one again obtains nicp niBu while in the pure envelope limit, 



mcP = '«env,P mp, 



(35) 



where mp = (pI^'^G^/^ / c^) M is the polytropic analog of of msi 
for a polytropic sphere. Note, once again, the factor of that 
enters since nicp is normalized to Ccore- 

The critical mass of the CPS is plotted in the isentropic case 
for various envelope indices and r, = 1 in Figure 6. Note the 
asymptotic approach to the Bonnor-Ebert mass of the core at 
ri/r, = 1, and the different limiting values (according to n) at 
'";/ rs = 0. The latter values, which are for pure envelopes, differ 
from those of Shu et al (1972) due to the difference in non- 
dimensionalization. Their values are normalized to the surface 
temperature so that the critical mass is always less than the 
Bonnor-Ebert mass (i.e. < Mbe = nisE for NIPs), whereas 
our values are normalized to the central temperature, and thus 
(for pure envelopes) mp^r can exceed msE- The critical mass of 
a pure polytrope has a maximum of mp„ ~ 1 .4 mpE at « = -2.90 
(7p = 0.655); this is in fact the maximum for all isentropic 
CPSs. CPSs with -2.12<n<-l have a maximum in mcp^a 
at intermediate r,/rs. In the limit of a constant pressure poly- 



trope, n 



-1 hp 



0), the critical mass shrinks to zero as 



the pressure gradient vanishes, while the mass for \n\ ^ 100 is 
essentially constant at the Bonnor-Ebert value for all ry/r,. 

In Figures la and b, we focus on the fiducial envelope index 
n = -3(7p = 2/3)in order to examine the behavior of the critical 
mass as a function of r,-. Note the similarity of these curves to 
those of the CIS, shown in Fig. 3. The "loops" appearing in Fig. 
la for r, > 3, which mark the end of the envelope's dominant 
mass contribution, indicate that configurations with the same 
fractional core radius can have different critical masses. As in 
the case of the CIS (§0), this multi-valuedness disappears in a 
plot of mcp^ci versus the fractional core mass, Mi/M, (Fig. lb). 

5.3. Velocity Dispersion 

The observable component of the sound speed, c(r), is the 
velocity dispersion, <j{r), given by equation (E2). The structure 
of the equilibrium CPS gives P{r) and p(r), so c(r) is easily cal- 
culated from equation (pl|). In order to gain a qualitative feel for 
a{r) and for the nonthermal velocity dispersion, a^nir) = a„{r), 
we set va = Q and plot these quantities for an isentropic CPS 
with n = -3, Ti= 1, and various values of rj/r, in Figure 8a. 
Also plotted are the same quantities for a pure n = -3 poly- 
trope, with the central dispersion normalized to Ccore- Outside 

the kinetic temperature may increase outward. Then, an assumed T^i^(r) gives (T,i,(r) 
1 as n — > —1, this is not true of (,Pc/Ps)ct', rather, (pc/ Ps)a — > 5.13 in this limit (MH). 
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of the thermal core, both dispersions increase monotonically 
with radius. Note how the slope of the nonthermal dispersion 
in the envelope increases with increasing core fraction. Fig. 8^? 
shows how, in certain cases (here for n = -1 .5, corresponding to 
7^ = 1 /3), both cr and its slope can exceed the pure polytropic 
values for a finite range of rz/rj (here 0.2 <J r,/rj <^ 0.55). This 
behavior will become increasingly important once we set out to 
model the observed increase of linewidth with size. 

In order to compare these model velocity dispersions with the 
observed molecular linewidths, one must first calculate the pro- 
jected linewidth through the region of the cloud having n > n^, 
where Wcr is the critical density for the excitation of a given 
molecular line (e.g. Genzel 1992). We will discuss this proce- 
dure more fully, and present detailed modehng, in a subsequent 
paper 

5.4. The Non-isentropic CPS 

The CPS offers an improvement over the CIS, in that it is 
qualitatively consistent with the observed increase of velocity 
linewidth with size in molecular clouds. But, as is evident from 
the critical masses plotted in Fig. 6, the isentropic CPS shares 
a common problem with the CIS in that it cannot quantitatively 
reproduce the observed mass range of PPCs; in order to obtain 
masses and density contrasts that are as large as those inferred, 
one requires r, ^ 3, a jump in for which there is neither the- 
oretical nor observational justification. As we have argued, in- 
sofar as molecular clouds can be described as polytropes, they 
are more likely to be non-isentropic in character It is therefore 
of some benefit to examine the quantitative properties of such 
polytropes. 

The theory of single-component, non-isentropic polytropes 
was reviewed and extended by MH. They showed that, for a 
given 7p (or «), the value of {pd Ps)cr increases as 7 increases 
above 7^, and reaches infinity at a well-defined value 7 = 700, 
given by 

327p(2->) 32(«2-l) 



7o 



(6->)' 



(5n-l)2 



(36) 



corresponding to the family of singular polytropic configura- 
tions. Observe that in the isothermal case (7^, = 1, n = ±cx)), 
7oo = 32/25, as found by Yabushita (1968). 

Again focussing on the fiducial value n = -3, in Figures 
9a and b we plot contours of maximum density contrast, 
iPc/ Pi)™"- for Ti = 1 and r,- = 2. Recall that these maxima oc- 
cur for different values of r,/rs, depending on the exact val- 
ues of T,, 7core, and 7env. The adiabatic index in the envelope, 
7env, ranges between the isentropic value, 7env = 7p,env = 2/3, 
and 7env = 7oo = 1, while 7core is confined to the correspond- 
ing range for the isothermal core (§4.3). Note the qualitative 
similarity between Fig. 9b and Figs. 4a and b. 

Since the CPS has one more free parameter than the CIS, 
it is interesting to examine properties of the critical state as n 
varies. In Figure lOa, we plot {pd Ps)f^^ and (Pf/f 5)™" versus 
7;j.env = ! + !/« for T; = 1 , 2, 3, and 5. Here, we fix 7coi-e = 1.1 and 
take 7env = 7;j.env + . 1 5 , in order to focus on the variation with n . 
Because 7env = 7p,env + 0.15 and equation ( |3^ have two points 
of intersection in the (7env,7;;.env) plane, there are two distinct 
values of 7p,env where (pc/Pi)™" iPc/Ps)^^^ become infi- 
nite: "fp.eny = 0.22 and 1.18. The 7p,env at which (pc/psY™^ is 
a minimum is roughly insensitive to r,; it lies in the the range 
0.75 ^ 7p,e„v ^ 0.90 (-10 ^ « ^ -4) for r, < 3. 

As in the case of the non-isentropic CIS, critical masses well 
in excess of the Bonnor-Ebert mass of the core may obtain for 



the non-isentropic CPS, even when r, = 1. This divergence of 
the mass is associated with a divergence of / pj'^\ if the prop- 
erties of the surface are known, then of course there is no di- 
vergence (§^. This behavior is illustrated in Fig. 10^?, which 
shows the maximum critical mass (once again, obtained from 
the run of nicp.a versus rz/r,) for CPSs of different r, , as a func- 
tion of 7y,.env (solid curvcs). As above, we fix 7core = 1.1. The 
dotted curve is for isentropic pure polytropes, for which mp„ 
is always finite. The dashed curves correspond to pure poly- 
tropes for the same values of r,, that is, having central tem- 
peratures = Ti Tcore. Notice how the CPS mass is always 
equal to or greater than the pure polytropic mass. Also, the 
former greatly exceeds the latter at low values of 7p.env, diverg- 
ing more quickly than the pure polytrope for all values of r,. 
Conversely, at large 7p.env the effect of the isothermal core is 
clearly seen, as ot™;^^. remains finite; by contrast, mp„ rapidly 
approaches zero as 7p.env ^6/5 and 7 exceeds 4/3 (MH). The 
limiting mass as 7p.env ^ 6/5 is slightly less than niBE simply 

because 7coie > 7/),core = 1 (§P^)- 

5.5. Mean Quantities in the CIS and CPS 

The existence of large pressure contrasts in the CIS and CPS 
does not necessarily imply that their mean pressures are also 
larger than the Bonnor-Ebert value. In fact, this is possible only 
for non-isentropic models. The mean pressure is given by 

P={d)-p, 

where 

Jo ' r^pdr 

is the mass-average of the squared sound speed, and p = 
M{r,)/{ATTrl/3) is the mean density. We have calculated the 
mean-to-surface pressure ratio, P/Ps, for a range of models, 
with the following results: (i) P/Ps of the critical isentropic 
CIS or CPS is always less than or equal to the Bonnor-Ebert 
value, (P/Ps)be = 2.43. (ii) The non-isentropic CIS has P/f, > 
(P/Ps) BE for all ri/rs, as shown in Figure 11a. If 7core > 7env, 
then the maximum P/Ps is obtained in the pure core limit, and 
vice-versa for 7env > 7core- Since the upper bound is for a single 
component isothermal sphere, we can use the results of MH to 
set an upper limit of P/P, < 3.78. (iii) The non-isentropic CPS 
has P/Ps less than the corresponding non-isentropic CIS. When 
the core size is significant (r,/ rj > 0.2, as shown in Fig. 1 1^7 for 
n = -3), the pressure ratio can substantially exceed {P/Ps)be\ 
otherwise, P/Ps ~ (P/Ps)nip < (P/P,)be- The last result fol- 
lows from the fact that P/Ps is not greatly enhanced for non- 
isentropic versus isentropic NIPs, even as 7 700. As shown 
by MH, 



< 1.26 



3(2- 



NIP 



(6-57;,) 



>^ = 1.26.^^"- 



1) 



(n-5) 



so that (P/Ps)max = 1.89 and 3.78 for n = -3 and n = -oo, re- 
spectively. These limits are consistent with the limiting values 
shown in Fig. lib. 

5.6. Summary — CPS 

The results of this section lead to two general conclusions 
concerning the approximation of mean and fluctuating mag- 
netic fields by isotropic pressure components (see also MH). 
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First, both the maximum dimensionless mass nicp that can be 
supported in a given model and the maximum pressure contrast 
are determined by the values of the polytropic and adiabatic in- 
dexes. As j-'jp (7 < 4/3) increases in either the core or the en- 
velope, {Pc/Ps)cp.a becomes orders of magnitude larger than its 
isentropic value. This behavior qualitatively mimics that found 
in equilibrium calculations of magnetized isothermal clouds, 
which feature stable pressure contrasts well in excess of the 
Bonnor-Ebert value for strong magnetic fields (Mouschovias 
1976; Tomisaka et al 1988). However, the behavior of the ratio 
of the mean pressure to that at the surface is quite different; we 
have extended the results of MH to composite polytropes with 
negative-index envelopes and shown that P/P, < 4, just as in 
the case of single component NIPs. 

Second, the restriction to isentropic NIPs, whether single- 
component or composite with r,- ~ 1, puts severe limitations 
on the mass range of the resulting models. Equilibrium clouds 
of mass substantially larger than the Bonnor-Ebert mass of the 
core do not exist under these restrictions. On the other hand, in 
the more likely case of non-isentropic polytropes (7env > Ip.em 
due to the presence of a mean magnetic field and Alfven waves), 
both single-component NIPs and composite CPSs can have 
much larger critical masses. 

The distinction between isentropic and non-isentropic poly- 
tropes is best seen by putting the critical mass in dimensional 
form. The mass of a composite polytrope is 



Mcp = , mcp = 2.2 Mr.^ 



2x lO^cm-^K 



1/2 



mcp, 



(37) 



Thus 2.6 Mq < Mrp < 3.6 Mq, assuming these representa- 
tive values (see eq ^5|). The larger upper limit, compared to the 
CIS (|34|), is due to the increased mass for pure n = -2.9 poly- 
tropes when expressed in terms of the central temperature. For 
a given surface temperature and pressure, Mcp.a < 3.6 Mq 
depends only weakly on whether the polytropic envelope is 
isentropic or not. However, the core temperature is rarely less 
than about 10 K; for a fixed core temperature and surface pres- 
sure, the maximum critical mass scales as cx (Pc/Pj)"^*'"'/'''''', 
which is limited for isentropic NIPs, but can be unbounded for 
non-isentropic ones. This fact has important consequences for 
the observed structure of embedded dense cores, a topic we take 
up in the following section. 

Although the composite polytropes discussed in this paper 
represent the simplest possible improvement on pure polytropic 
models, they already suggest a number of interesting implica- 
tions for star formation generally. We discuss a few of these 
below. 

6. DISCUSSION 

6.1. Near-Critical Cores with Massive Envelopes 

An interesting consequence of our results for the non- 
isentropic CPS (§5.4) is the existence of equilibria featuring 



The core, envelope, and total dimensionless masses of a 
n = -3 CPS with 7core = 1-2, 7env = 0.94 are displayed as a func- 
tion of the core fraction rj /r, in Figure 12fl. The critical mass of 
the core alone, mcme.cr, has also been plotted, in the same non- 
dimensionalization. For a fixed central sound speed Qoie and 
outer surface pressure Pj, the pressure on the core, P,, increases 

-1/2 

as the core fraction shrinks; hence, mcore.cr Mcore. cr « P, 
decreases with decreasing rz/r,. Fig. 12fl has a very interesting 
property: the mass of the core remains very near (and indeed, 
somewhat exceeds^) its critical mass for r, /r, > 1.3 x 10"^ de- 
spite the fact that the corresponding mass of the envelope varies 
by more than an order of magnitude over the same range. If 
observed in a high-density tracer, such a core would appear 
to have a virial parameter of nearly the Bonnor-Ebert value 
(a = 2.06; Bertoldi & McKee 1992), as is observed. In terms of 
our dimensionless variables, the virial parameter becomes (see 



1/2 



a = — = 47r — 

GM \ Ps 



Taking (cr^) = ^ = and m = mcore, we find acoi-e — 2 for 
all rz/r, > 1.3 x 10"^, as shown by the upper curve in Fig. 
\2b. It is of interest to compare Ocore with the virial parame- 
ter for the entire cloud, We estimated the mean velocity 
dispersion of the latter using equation (^2|), and by assuming 
(cr^) 3> (T(h (valid for r,- ^ r^) and an Alfven Mach number of 1, 



I.e. v'1 



:3(-2\ 



massive envelopes surrounding approximately Bonnor-Ebert 
mass cores, a property reminiscent of the observations reviewed 
in §1.2. Since the derived mass ratio of '^^CO envelopes to NH3 
cores is of order 10-100, we choose an example CPS which 
illustrates this mass range. 

^Recall that m > ma for non-isentropic states with (Pc/Ps)be < Pc/Pi < {Pc/Pi)t 



-^[a") (Myers & Goodman 1988; Crutcher 1999). This 
gives (cr^) ~ c^/3. Strictly speaking, the latter should change 
with Tj/rj, but this estimate suffices to get a rough idea of the 
magnitude of oftot- Finally, taking ^ = and m = mcp, we find 
that Qftot varies between 0.73 and 1.5 as ri/rg ranges from to 1 
(Fig. 12b). Values of cttot < 1 simply reflect the limited validity 
of our estimated (cr^) — as r, r,, the cloud becomes thermally 
dominated with (a^) — > c^, not c^/3 as we have assumed. 

Insofar as the derived virial parameters of many clumps 
on sub-GMC scales are much greater than unity (§1.2), we 
conclude (with Bertoldi & McKee 1992), that such clumps 
are pressure-confined. Sub-critical CPSs having M^w =10 — 
100 Mcore are easily constructed. In our view, these results 
present a reasonable explanation for why it is that PPC s, which 
are nearly critical by virial theorem estimates (§1.2), can ex- 
ist within clumps and GMCs that are far more massive but 
nonetheless stable or close to critical. This constitutes an ex- 
plicit realization of the formalism developed by Bonazzola et 
al (1987), who proposed that a turbulent cloud could be sta- 
ble on large scales and unstable on small ones (note that our 
models are stable by construction, but would become unstable 
if the core mass were increased further for rz/r, > 1.3 x 10"^). 
A more detailed comparison of this aspect of the models with 
existing data is, however, beyond the scope of this paper. 

6.2. The Formation of Very Low-Mass Stars 

The enhanced pressure contrasts attainable in the composite 
models may also have implications for forming very low-mass 
stars such as brown dwarfs. A simple argument for an isolated 
core forming a single star serves to illustrate this point. In or- 
der to form a star of mass M*, it is necessary to collect together 
at least a Bonnor-Ebert mass, Mbe = 2.6 (r/10 Kf {P,/{2 x 
10"* ks) K cm"^]"'/2 Mq, of gas. Observed cores in GMCs 

: (Fig. 1). 
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have surface pressures ranging from Ps/ks ~ 4 x 10** (Tau- 
rus, Cepheus) to 3 x 10^ K cm^ (Orion) (Solomon et al 1987; 
Bertoldi & McKee 1992), so it is easy to form stars with 
M» ~ 0.7- 1.8 Mq. But to create a star of much smaller mass 
(brown dwarfs have <^ 0.08Mq), it is necessary to either 
reduce the temperature below 10 K (which is difficult due to 
radiative trapping- de Jong et al 1980), or to increase Ps/ks 
well beyond the observed range. Inserting Mj < 0.08 Mq into 

the above, one finds Ps/kg > 1.3 x 10*^ K cm-*, nearly three 
orders of magnitude above the observed pressures. In a CPS 
with ffJcore = 0.08 Mq, howcvcr, it is the pressure on the core- 
envelope interface, P,, that counts. Here one is not limited by 
the Bonnor-Ebert value Pi/P,, = 14.04, since the pressure drop 
through a non-isentropic envelope can be much larger (Fig. 
10a). Thus, non-isentropic CPSs lead to the possibility of a 
structure near the brink of instability with a core small enough 
to make a brown dwarf. Whether it is possible for the core to 
undergo collapse without having a significant mass from the 
envelope fall on it, thereby increasing its mass, cannot be deter- 
mined from hydrostatic calculations such as ours. 

6.3. Protostellar Collapse 

One of the more robust predictions of studies of cloud col- 
lapse is the mass accretion rate onto the nascent central object. 
The well-studied SIS, for example, features a time-independent 
accretion rate, Msis = 0.975 cTj^G (Shu 1977). Observational 
estimates of M have been inferred for many Class and Class I 
infrared sources, from the peak of their spectral energy distribu- 
tions in the far-infrared (Hartmann 1998). These fall within the 
range 1.6 x IQ-^Mq yr"' (Class I) ;$M ^ IQ-'^Mq yr"' (Class 

0). For a typical fith ~ 0.2 km s~', the SIS is in accord with the 
lower limit of this range, but is inconsistent with the high M's 
of the Class objects. This again raises the question of the rele- 
vance of the SIS for the earliest stages of protostellar evolution 
(Foster & ChevaHer 1993; Safier et al 1997). 

We now attempt to predict how the collapse of a compos- 
ite polytrope will differ from that of previous models. Since 
we have modeled the core of a CPS with an isothermal sphere, 
the early behavior is likely to be similar to that of the Bonnor- 
Ebert sphere, which was investigated most recently by Foster 
& Chevalier (1993). They found that M is largest at the onset 
of collapse, then decreases, eventually approaching Ms/5. The 
late-time behavior is more uncertain, and depends on the char- 
acteristics of the envelope. If the envelope can be described by a 
NIP, then two recent studies suggest that the accretion rate will 
increase at late times: McLaughlin & Pudritz (1997) consid- 
ered expansion wave solutions in singular polytropic spheres 
and found that M cx t~^^" (n < -1); i.e., the accretion rate is 
small at early times, then increases mono tonic ally. Henriksen, 
Andre, & Bontemps (1997) studied the collapse of a pressure- 
less model (so that there is no expansion wave) with a piecewise 
continuous density profile: p oc r" in the core and p cx ^-^"/("-i) 
(corresponding to a singular polytropic sphere) in the envelope. 
Their results show that after an initial burst of accretion (M is 
formally infinite in their model, due to the pressureless, flat in- 
ner region), models with « < -1 pass through a minimum in 
M(> Msis)^ before increasing monotonically at large t as t~^l", 
in agreement with the results of McLaughlin & Pudritz. On the 
other hand, if the envelope is supported by a static magnetic 
field so that the late-time accretion is regulated by ambipolar 
diffusion, the accretion rate is similar to that of an SIS (Safier 
et al 1997). More definitive results for the collapse of a CPS 



await detailed numerical calculations. 

7. SUMMARY 

In this paper, we have developed a polytropic theory of 
molecular clouds that improves on previous work in several 
ways. The consideration of spatially separate core and envelope 
regions is in better accord with the known velocity structure of 
star-forming cloud cores. Such an approach allows a thermally- 
dominated core to coexist with a nonthermally-dominated en- 
velope in a single, self-gravitating equilibrium structure. In 
the formalism developed here, the core and the envelope are 
each represented by a single pressure component with prop- 
erties chosen to approximate the combined effects of thermal 
pressure and static magnetic fields in the core (assumed to be 
described by a polytropic index 7p =1), and these plus turbu- 
lent motions in the envelope. The results for isothermal en- 



velopes are summarized in §4.4, and those for polytropic en 



velopes in §5.6. Just as composite polytropic models of stars 



were created to reflect changing physical conditions with ra- 
dial scale, cloud models can now do the same. Our models are 
non-singular, a feature in better concord with observations of 
the densest pre-protostellar cores than the singular hypothesis. 
This agreement carries over to the quantitative comparison of 
both the density and velocity dispersion profiles with those ob- 
served, as we present in detail in a subsequent paper. 

The increase of temperature near the surface and of line 
width with size had led previous workers to consider negative- 
index polytropes (NIPs) as models of interstellar clouds (Shu et 
al 1972; Maloney 1988). In general, these models were based 
on the assumption that the gas is isentropic — i.e., that the poly- 
tropic index 7^, is the same as the adiabatic index 7. The max- 
imum pressure and density drops in such NIPs are always less 
than that of a Bonnor-Ebert sphere, however, suggesting that 
they are inadequate as models for molecular clouds. Isentropic 
composite polytropes of reasonable core fractions already pro- 
vide pressure and density enhancements well in excess of those 
of the isothermal sphere. The critical mass and radius of an 
isentropic composite polytropic sphere (CPS) can also exceed 
the Bonnor-Ebert mass of the core in cases where the kinetic 
temperature of the envelope is higher than that of the core. 

However, the pressure components in molecular gas, partic- 
ularly magnetic fields and turbulence (which we model with 
Alfven waves), are not isentropic (MH). Non-isentropic CPSs 
permit pressure and density drops in stable structures that are 
as large as those observed, and in addition permit masses that 
can be large compared with the Bonnor-Ebert mass of the core, 
even without a temperature difference between core and enve- 
lope. However, as found by MH for individual NIPs, the ratio 
of the mean pressure to the surface pressure for composite poly- 
tropes with negative-index envelopes must be less than a factor 
~ 4, which appears to be a general feature of hydrostatic models 
for self-gravitating clouds with outwardly increasing tempera- 
tures. The greater flexibility of the CPS allows one to transcend 
the limitations of many previous polytropic models, which were 
unable to reproduce even the bulk physical properties of PPCs, 
not to mention the detailed spatial structure of these quantities. 
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APPENDIX 

ON THE STABILITY CRITERION FOR POLYTROPIC 
MODELS OF MOLECULAR CLOUDS 



As discussed in § 1 .4 , the stability analysis of bounded poly- 



tropes has been pursued in the literature according to two very 
different assumptions about the response of the cloud to small 
perturbations. Some authors have assumed that the central tem- 
perature remains constant during the perturbation, whereas we 
(in agreement with MH) have assumed that the entropy remains 
constant. For an isothermal gas, the two assumptions are equiv- 
alent; but in general they are not. The issue is not one of mere 
formalism: whether a given model is stable or not is highly sen- 
sitive to the choice. 

The constant central temperature condition was justified by 
Maloney (1988) on the basis of the ". . .general uniformity of 
gas temperatures observed in molecular clouds over a large 
range of spatial and density scales ..." As Turner et al (1992) 
pointed out, however, it is the value of for a given cloud dur- 
ing a perturbation that is relevant to the issue of stability. Still 
asserting the constancy of Tc, however, the latter authors argued 
that it ". . . is determined by the cloud's extinction and the local 
UV [field], and variations in [71] due to (small) variations in 
[the boundary pressure] may reasonably be considered negligi- 
ble." 

The issue becomes more involved when we consider that 
these authors view the increasing outward temperature of NIPs 
as being partly or entirely nonthermal in origin. In the models 
of both Maloney and McLaughlin & Pudritz (1996), e.g., the 
kinetic temperature is spatially constant (and supposed equal to 
Tc), but superposed on this is an isotropic, turbulent velocity 
field (5v, whose effective temperature, Teff oc {S\)^, behaves like 
that of a NIP. This is, in fact, in accord with recent observations 
of Bok globules showing that the gas temperature throughout 
a given globule is remarkably uniform, even though the non- 
thermal linewidth can be up to twice the thermal value (Lemme 
et al 1996). The origin of the turbulent motions producing the 
nonthermal linewidths is the object of much speculation (e.g. 
Alfven waves, cloud-cloud motions, protostellar outflows, etc.), 
but it is certainly different than the mechanisms that determine 
the thermal motions in the gas. 

In a polytropic model, the surface temperature is related to 
the central temperature by 



nonthermal velocity dispersion anywhere in the cloud. Further- 
more, as pointed out by MH, both Maloney and McLaughlin 
& Pudritz considered isentropic polytropes. Such polytropes 
cool upon compression if 7^ < 1, as they assumed. In order 
to maintain a constant central temperature, energy would have 
to be supplied, and no plausible energy source has been identi- 
fied. In the language of §2.1 and MH, the isothermal perturba- 
tions considered by these authors are neither locally nor glob- 
ally adiabatic, since heat is not simply redistributed but rather 
is injected by some unknown agent. 

In our analysis, we avoid these contradictions and follow 
most previ ous a uthors by taking the perturbations to be locally 
adiabatic (§2.1): no heat enters or leaves the system (cf. Shu 
et al. 1972; Viala & Horedt 1974; Stabler 1983). However, 
we depart from these treatments by allowing for non-isentropic 
pressure components, whose adiabatic index 7, which dictates 
the stability properties of the gas, differs from the polytropic 
index •jp [or n = 1/(1 -7^)], which determines the equilibrium 
structure. For NIPs in particular, taking 7 > jj, moves the 
critical point of equilibrium to larger dimensionless radii, and 
more centrally condensed states (Fig. 1), than in the isentropic 
(7 = jp) case. These newly-accessible equilibria are identical to 
those assumed to exist by the above authors, but are not derived 
from the questionable assumption of unnamed energy sources. 



Tr 



(38) 



The central region is often sufficiently shielded that the temper- 
ature there is determined by cosmic ray heating. Thus, while 
it is plausible that Tc remains constant during the perturbation, 
there is no justification for assuming that 7^, which is repre- 
sented in the turbulent models by the velocity dispersion at the 
surface of the cloud, will change just so as to satisfy this relation 
with rc=const. Indeed, there is no known relation between the 
physical processes determining the kinetic temperature and the 
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log 10 (Pc/Ps) 

Fig. 1. — Equilibrium sequences of the isotliermal sphere (curves with filled circles) and an n = -3 polytrope (curves with open circles). Solid lines show the 
dimensionless mass; dashed lines, the dimensionless radius. Dots indicate critical points at particular values of 7, as labeled. 
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Fig. 2. — Density and pressure contrasts of the critically stable, isentropic CIS. (a) Critical density contrast as a function of fractional core radius, n/rs. Each 
point on a curve represents a single, critically stable, equiUbrium model. Solid curves are center-to-surface contrast, dashed curves center-to-core boundary contrast. 
Curves are labeled by values of t;. The filled circle on the ordinate indicates a density contrast of Ti(pc/ps)\g = 50(14.04)^ = 9855. (b) Critical density contrast as 
a function of fractional core mass, using same notation as in (a), (c) Critical pressure contrast as a function of fractional core radius, using same notation as in (a). 
The filled circle on the ordinate indicates a pressure contrast of (Pc/PsigE = (14.04)^ = 197. 
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Fig. 4. — Maximum pressure contrast of the critically stable, non-isentropic CIS. Contours of equal logio(Pc/^j)S'" ^re plotted as a function of the core and 
envelope adiabatic indices, 7core and 7env, respectively, for (a) r, = 2; (b) t; = 3. 
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Fig. 5. — Pressure and density contrasts of the critically stable, isentropic CPS, in the cases (a) t, = 1; (b) t, = 2; and (c) r,- = 3. Solid curves show the center-to- 
surface density contrast; dashed curves, the center-to-surface pressure contrast. Both are plotted as a function of the fractional core radius. Curves are labeled by 
values of n. 
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Fig. 7. — Dimensionless mass of the critically stable, isentropic CPS. (a) Critical mass, normalized to the Bonnor-Ebert mass mg£, as a function of fractional core 
radius, r,/ r,. Solid curves show the core mass alone; dashed curves, the envelope mass; and heavy solid curves the sum of the two. Curves are labeled by values of 
Tj. (b) Critical mass as a function of fractional core mass, using same notation as in (a). 
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Fig. 8. — Velocity dispersion in the critically stable, isentropic CPS, assuming va = 0. (a) Total (top curves) and nonthermal (bottom curves) velocity dispersions 
for CPSs with r, = 1, n = -3, and various values of r,/ri. From left to right: rj/r., = 0.2, 0.4. 0.6, and 0.8. Also shown are the corresponding dispersion for a pure 
« = -3 polytrope (dashed curves), (b) Similar to (a), but for n = -1.5. Only the total dispersion is plotted. 




Fig. 9. — Maximum density contrast of the critically stable, n = -3 non-isentropic CPS. Contours of equal logio(Pc//'s)S™ ^re plotted as a function of the core 
and envelope adiabatic indices, 7core and 7env, respectively, for (a) t, = 1; (b) t, - 2. 
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7p.a„v = 1 + 1/n 7p,.„ = 1 + 1/n 

Fig. 10. — Properties of non-isentropic CPSs with 7core = 1.1 and ■y^m = 7p.cnv + 0.15. (a) Maximum critical pressure (dashed curves) and density (solid curves) 
contrasts as a function of 7p,eiiv. Curves are labeled by r, . (b) Maximum critical masses as a function of 7;),eiiv (solid curves). Also shown are the critical masses of 
pure polytropes with Tcore = t,- rcore(T,- = 1) (dashed curves). The dotted line shows the mass of isentropic polytropes. 




Fig. 11. — (a) Mean pressure contrasts (soUd curves) as a function of fractional core radius in the t; = 1, non-isentropic CIS. Each curve is labeled by its 
corresponding value of (7core , 7eiiv). The dotted line shows the mean sound speed contrast, which is equal to unity since r; = 1. (b) Mean pressure contrasts (soUd 
curves) in the t,- = 1, non-isentropic CPS with n = -3 and 7env = 0.85. Mean density contrasts (dashed curves) and mean sound speed contrasts (dotted Unes) are also 
shown. Within each set of curves, the corresponding value of 7core is, from bottom to top, 1, 1.1, and 1.2. 
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Fig. 12. — (a) Dimen.sionless masses of a critically stable, non-isentropic CPS with r, = 1, as a function of the fractional core radius. The solid curve shows the 
core mass alone; the dashed curve, the envelope mass; and the heavy solid curve the total of the two. The dotted curve is the critical mass of the core alone, (b) 
Virial parameters for the core alone (solid curve) and the entire cloud (heavy solid curve). See text for details. 



